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Abstract. We study the noise of the cotunneling current through one or several tunnel- 
coupled quantum dots in the Coulomb blockade regime. The various regimes of weak and 
strong, elastic and inelastic cotunneling are analyzed for quantum-dot systems (QDS) 
with few-level, nearly-degenerate, and continuous electronic spectra. In the case of weak 
cotunneling we prove a non-equilibrium fluctuation-dissipation theorem which leads to 
a universal expression for the noise-to-current ratio (Fano factor). The noise of strong 
inelastic cotunneling can be super-Poissonian due to switching between QDS states car- 
rying currents of different strengths. The transport through a double-dot (DD) system 
shows an Aharonov-Bohm effect both in noise and current. In the case of cotunneling 
through a QDS with a continuous energy spectrum the Fano factor is very close to one. 



1. Introduction 

In recent years, there has been great interest in the shot noise in mesoscopic 
systems [1], because it contains additional information about correlations, 
which is not contained, e.g., in the linear response conductance. The shot 
noise is characterized by the Fano factor F = S/el, the dimensionless 
ratio of the zero-frequency noise power S to the average current /. While 
it assumes the Poissonian value F = 1 in the absence of correlations, it 
becomes suppressed or enhanced when correlations set in as e.g. imposed 
by the Pauli principle or due to interaction effects. In the present paper 
we study the shot noise of the cotunneling [2, 3] current. We consider the 
transport through a quantum-dot system (QDS) in the Coulomb blockade 
(CB) regime, in which the quantization of charge on the QDS leads to a 
suppression of the sequential tunneling current except under certain res- 
onant conditions. We consider the transport away from these resonances 
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and study the next-order contribution to the current 1 (see Fig. 1). We 
find that in the weak cotunneling regime, i.e. when the cotunneling rate 
I/e is small compared to the intrinsic relaxation rate w- m of the QDS to its 
equilibrium state due to the coupling to the environment, I/e> w in , the 
zero-frequency noise takes on its Poissonian value, as first obtained for a 
special case in [6]. This result is generalized here, and we find a universal 
relation between noise and current for the QDS in the first nonvanishing 
order in the tunneling perturbation. Because of the universal character of 
this result Eq. (12) we call it the nonequilibrium fluctuation-dissipation 
theorem (FDT) [7] in analogy with linear response theory. 

One might expect however that the cotunneling, being a two-particle 
process, may lead to strong correlations in the shot noise and to the de- 
viation of the Fano factor from its Poissonian value F = 1. We show in 
Sec. 4 that this is indeed the case for the regime of strong cotunneling, 
//e» W{ n . Specifically, for a two-level QDS we predict giant (divergent) 
super- Poissonian noise [8] (see Sec. 5): The QDS goes into an unstable mode 
where it switches between states 1 and 2 with (generally) different currents. 
In Sec. 6 we consider the transport through a double-dot (DD) system as an 
example to illustrate this effect (see Eq. (37) and Fig. 2). The Fano factor 
turns out to be a periodic function of the magnetic flux through the DD 
leading to an Aharonov-Bohm effect in the noise [9]. In the case of weak 
cotunneling we concentrate on the average current through the DD and find 
that it shows Aharonov-Bohm oscillations, which are a two-particle effect 
sensitive to spin entanglement. 

Finally, in Sec. 7 we discuss the cotunneling through large QDS with 
a continuum spectrum. In this case the correlations in the cotunneling 
current described above do not play an essential role. In the regime of 
low bias, elastic cotunneling dominates transport, [2] and thus the noise is 
Poissonian. In the opposite case of large bias, the transport is governed 
by inelastic cotunneling, and in Sec. 7 we study heating effects which are 
relevant in this regime. 

2. Model system 

In general, the QDS can contain several dots, which can be coupled by tun- 
nel junctions, the DD being a particular example [6]. The QDS is assumed 
to be weakly coupled to external metallic leads which are kept at equilib- 

1 The majority of papers on the noise of quantum dots consider the sequential tun- 
neling regime, where a classical description ( "orthodox" theory) is applicable [4] . In this 
regime the noise is generally suppressed below its full Poissonian value F = 1. This 
suppression can be interpreted [5] as being a result of the natural correlations imposed 
by charge conservation. 
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rium with their associated reservoirs at the chemical potentials m, I = 1,2, 
where the currents l\ can be measured and the average current / through 
the QDS is defined by Eq. (5). Using a standard tunneling Hamiltonian 
approach [10], we write 



H = H + V, H = H L + H s + H int , (1) 

h l= ^2^2 £ k4k c ik , h s = ^2 £ P d l d p > ( 2 ) 

1=1,2 k V 

V = Y, (A + D\), A = £ T lkp c\ k d p , (3) 

1=1,2 k,p 



where the terms Hl and H$ describe the leads and QDS, respectively (with 
k and p from a complete set of quantum numbers), and tunneling between 
leads and QDS is described by the perturbation V. The interaction term 
H mt does not need to be specified for our proof of the universality of noise 
in Sec. 3. The iV-electron QDS is in the cotunneling regime where there is 
a finite energy cost A±(l, N) > for the electron tunneling from the Fermi 
level of the lead I to the QDS (+) and vice versa (— ). This energy cost is of 
the order of the charging energy Ec and much larger than the temperature, 
A±(l,N) ~ Eq S> ksT, so that only processes of second order in V are 
allowed. 

Tj T 2 



A+(1.N) 



Hi 



lead 1 

Figure 1. The quantum dot system (QDS) is coupled to two external leads I = 1, 2 via 
tunneling barriers. The tunneling between the QDS and the leads is parametrized by the 
tunneling amplitudes T;, where the lead and QDS quantum numbers k and p have been 
dropped for simplicity, see Eq. (3). The leads are at the chemical potentials jJ.1,2, with an 
applied bias A/i = hi — [ii. The eigenstates of the QDS with one added electron (N + 1 
electrons in total) are indicated by their energies E\, E2, etc., with average level-spacing 
8E. In the cotunneling regime there is a finite energy cost A±(l, N) > for the electron 
tunneling from the Fermi level of the lead I to the QDS (+) and vice versa (— ), so that 
only processes of second order in V (visualized by two arrows) are allowed. 
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To describe the transport through the QDS we apply standard meth- 
ods [10] and adiabatically switch on the perturbation V in the distant 
past, t = to — ► — oo. The perturbed state of the system is described 
by the time-dependent density matrix p(t) = e~ iH ( t ~ t °^ poe iH ( t ~ t °\ with 
po being the grand canonical density matrix of the unperturbed system, 
po = Z~ 1 e~ K / kBT , where we set K = Hq — ^ mN\. Because of tunneling 
the total number of electrons in each lead Ni = Y^k c lk°lk ^ s 110 l° n g er 
conserved. For the outgoing currents Ii = eiVj we have 

I l = ei[V,N l ] = ei{D\-Di). (4) 

The observables of interest are the average current I = I2 = —I\ through 
the QDS, and the spectral density of the noise Sw(uj) = f dtSw(t) exp(icut), 

J, = Trp(0)4 S w (t) = ReTrp(0)^(t)<5^(0) , (5) 

where 5Ii — Below we will use the interaction representation where 
Eq. (5) can be rewritten by replacing p(0) — > po and Ii(t) — > U^(t)Ii(t)U(t), 
with 



U(t) = Texp 



-i I dt'V(t') 

J —CO 



(6) 



In this representation, the time dependence of all operators is governed by 
the unperturbed Hamiltonian Hq. 



3. Weak cotunneling: Non-equilibrium fluctuation-dissipation 
theorem 



In this section we prove the universality of noise of tunnel junctions in 
the weak cotunneling regime I/e <C w- m keeping the first nonvanishing 
order in the tunneling Hamiltonian V. Since our final result (12) can be 
applied to quite general systems out-of-equilibrium we call this result the 
non-equilibrium fluctuation-dissipation theorem (FDT). In particular, the 
geometry of the QDS and the interaction H ult are completely arbitrary for 
the discussion of the non-equilibrium FDT in this section. 

We note that the two currents are not independent, because I2] 7^ 
0, and thus all correlators Su> are nontrivial. The charge accumulation 
on the QDS for a time of order A^ 1 leads to an additional contribution 
to the noise at finite frequency u. Thus, we expect that for uj ~ A± the 
correlators Sw cannot be expressed through the steady-state current I only 
and thus / has to be complemented by some other dissipative counterparts, 
such as differential conductances Gw. On the other hand, at low enough 
frequency, to <C A±, the charge conservation on the QDS requires 5I S = 
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(512 + 5I\)/2 « 0. Below we concentrate on the limit of low frequency and 
neglect contributions of order of w/A± to the noise power. In the Appendix 
we prove that S ss ~ (uj/ A±) 2 (see Eq. (65)), and this allows us to redefine 
the current and the noise power as I = Id = (l2—Ii)/2 and S(oj) = Sddi 1 ^)- 2 
In addition we require that the QDS is in the cotunneling regime, i.e. the 
temperature is low enough, ksT <C A±, although the bias A/j, is arbitrary 
as soon as the sequential tunneling to the dot is forbidden, A± > 0. In 
this limit the current through a QDS arises due to the direct hopping 
of an electron from one lead to another (through a virtual state on the 
dot) with an amplitude which depends on the energy cost A± of a virtual 
state. Although this process can change the state of the QDS (inelastic 
cotunneling), the fast energy relaxation in the weak cotunneling regime, 
Wm 3> I/e, immediately returns it to the equilibrium state (for the opposite 
case, see Sec. 4). This allows us to apply a perturbation expansion with 
respect to tunneling V and to keep only first nonvanishing contributions, 
which we do next. 

It is convenient to introduce the notation Di{t) = dt' Di(t'). We 
notice that all relevant matrix elements, (N\Di(t)\N + 1) ~ e~ lA+t , (N — 
l\Di(t)\N) ~ e tA - t , are fast oscillating functions of time. Thus, under 
the above conditions we can write = 0, and even more general, 

I-00 Di(t)e ±tujt = (note that we have assumed earlier that uj <C A±). 
Using these equalities and the cyclic property of the trace we obtain the 
following results (for details of the derivation, see Appendix A), 

00 

I = e J dt([A*(t),A(0)]), A = D 2 d\ + d\d 2 , (7) 

—00 

00 

S(lo) = e 2 J dt cos(^)({^ t (t),^l(0)}) , (8) 

—00 

where we have dropped a small contribution of order u/A± and used the 
notation (...) = Trpo(- • ■)• 

Next we apply the spectral decomposition to the correlators Eqs. (7) 
and (8), a similar procedure to that which also leads to the equilibrium 
fluctuation-dissipation theorem. The crucial observation is that [Hq,Ni] = 
0, I = 1, 2. Therefore, we are allowed to use for our spectral decomposition 

2 We note that charge fluctuations, SQ(t) = 2 dt'5I s (t'), on a QDS are also relevant 
for device applications such as SET [11]. While we focus on current fluctuations in the 
present paper, we mention here that in the cotunneling regime the noise power {SQ 2 }^ 
does not vanish at zero frequency, (5<5 2 ) tJ= o = 4cj _2 5 ss (o;)|„^o 7^ 0. Our formalism is 
also suitable for studying such charge fluctuations; this will be addressed elsewhere. 
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the basis |n) = \E n , N\, N%} of eigenstates of the operator K = Hq — 
^2i/J>iNi, which also diagonalizes the grand-canonical density matrix po, 
p n = (n|po|n) = Z^ 1 exp[— E n /k B T]. We introduce the spectral function, 



A(u) = 2ir^2(p n + p m )|(m|^|n)| 2 <5(a; + E n - E m ) , 



(9) 



and rewrite Eqs. (7) and (8) in the matrix form in the basis |n) taking 
into account that the operator A, which plays the role of the effective 
cotunneling amplitude, creates (annihilates) an electron in the lead 2 (1) 
(see Eqs. (3) and (7)). We obtain following expressions 



I(Ap) = etanh 



Ap 
2k B T 



A(Ap), 



S(u;,Ap) = ^-^2A(Ap±u J ) 



(10) 
(11) 



± 



We note that because of additional integration over time t in the amplitude 
A (see Eq. (7)), the spectral density A depends on p\ and p2 separately. 
However, away from the resonances, uj <C A±, only A^-dependence is 
essential, and thus A can be regarded as being one-parameter function. 3 
Comparing Eqs. (10) and (11), we obtain 



S(w, ^) = ^Y1 coth 



± 



Ap ± uj 



2k B T 



I(Ap±uj) 



(12) 



up to small terms on the order of u/A±. This equation represents our 
nonequilibrium FDT for the transport through a QDS in the weak cotun- 
neling regime. A special case with T, oj = 0, giving S = el, has been derived 
earlier [6] . To conclude this section we would like to list again the conditions 
used in the derivation. The universality of noise to current relation Eq. (12) 
proven here is valid in the regime in which it is sufficient to keep the first 
nonvanishing order in the tunneling V which contributes to transport and 
noise. This means that the QDS is in the weak cotunneling regime with 
u, k B T < A±, and I/e <C w in . 



4. Strong cotunneling: Correlation correction to noise 

In this section we consider the QDS in the strong cotunneling regime, W[ n 
I/e. Under this assumption the intrinsic relaxation in the QDS is very slow 

3 To be more precise, we neglect small w-shift of the energy denominators A±, which 
is equivalent to neglecting small terms of order u>/A± in Eq. (11). 
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and will in fact be neglected. Thermal equilibration can only take place via 
coupling to the leads (see Sec. 7). Due to this slow relaxation in the QDS we 
find that there are non-Poissonian correlations AS in the current through 
the QDS because the QDS has a "memory" ; the state of the QDS after the 
transmission of one electron influences the transmission of the next electron. 
The microscopic theory of strong cotunneling has been developed in Ref. [5] 
based on the density-operator formalism and using the projection operator 
technique. Here we discuss the assumptions and present the results of the 
theory, equations (14), (15), and (17-19), which are the basis for our further 
analysis in the Sees. 5 and 6. 

First, we assume that the system and bath are coupled only weakly 
and only via the perturbation V, Eq. (3). The interaction part H mt of the 
unperturbed Hamiltonian Hq, Eq. (1), must therefore be separable into 
a QDS and a lead part, H[ nt = H 1 ^ + H™ 1 . Moreover, Ho conserves the 

number of electrons in the leads, [Hq,Ni] = 0, where Ni = Ylk c ik Clk ~ r ^^ ie 
assumption of weak coupling allows us to keep only the second-order in V 
contributions to the "golden rule" rates (15) for the Master equation (14). 

Second, we assume that in the distant past, to — ► — oo, the system is in 
an equilibrium state 

Po = Ps®p L , PL = ^e~ K ^ kBT , (13) 

Zh 

where Zl = Tr exp[— K^/ksT], Kl = -Hz, — X)z piNi, and pi is the chemical 
potential of lead I. Note that both leads are kept at the same temperature 
T. Physically, the product form of po in Eq. (13) describes the absence 
of correlations between the QDS and the leads in the initial state at to- 
Furthermore, we assume that the initial state po is diagonal in the eigenbasis 
of Hq, i.e. that the initial state is an incoherent mixture of eigenstates of 
the free Hamiltonian. 

Finally, we consider the low-frequency noise, uj <C A±, i.e. we neglect 
the accumulation of the charge on the QDS (in the same way as in the 
Sec. 3). Thus we can write Su(oj) = —Si^i^lo) = S(u>). This restriction will 
be lifted in the end of the Sec. 6.1. 

We note that the above assumptions limit the generality of the results 
of present section as compared to those of Sec. 3. On the other hand, they 
allow us to reduce the problem of the noise calculations to the solution of 
the Master equation 

Pn(t) = ^ [W nm p m {t) - W mn p n (t)] , (14) 

m 

with the stationary state condition {wnmPm ~ w mn p n ) = 0. This "clas- 
sical" master equation describes the dynamics of the QDS, i.e. it describes 
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the rates with which the probabilities p n for the QDS being in state \n) 
change. The rates w nm = Ylii'=i 2 w nm{V , are the sums of second-order 
"golden rule" rates 

w nm (l\l)=2n^\(n\(DlD l ,)\ m )\H(E m -E n -Afi U i)pL, m . (15) 

m,n 

for all possible cotunneling transitions from lead I to lead I'. In the last 
expression, Apw = pi — pi> denotes the chemical potential drop between 
lead / and lead I', and pL,m = {m\pL\m). We have defined the second order 
hopping operator 

(D},D 1 ,) = D 1 ,D} + d}d v , (16) 

where Di is given in Eq. (3), and Di = f ^ Di(t)dt. Note, that (Dj, Dy) is 
the amplitude of cotunneling from the lead I to the lead /' (in particular, 
we can write A = (D\,D2), see Eq. (7)). The combined index m = (m,fh) 
contains both the QDS index m and the lead index fh. Correspondingly, 
the basis states used above are |m) = \m)\m) with energy E m = E m + E m , 
where \m) is an eigenstate of H$ + H 1 ^ with energy E m , and \m) is an 
eigenstate of Hl + Hf^ — piNi with energy E m . 

For the average current / and the noise power S(u>) we obtain [5] 

I = eJ2 W irnPm, = W nm (2, 1) - W nm (l, 2), (17) 

mn 

S(lo) = e 2 ^K m (2 1 l)+ Wnm (l 1 2)]p m + A%), (18) 

mn 

AS(uj) = e 2 J2 wI nm S Pmn>(uWn> m >Pm>, (19) 
n,m,n',m' 

where 8p nm (uj) = p nm (uj) — 2irS(uj)p n , and p n is the stationary density 
matrix. Here, p nm {uS) is the Fourier-transformed conditional density matrix, 
which is obtained from the symmetrized solution p n (t) = p n (—t) of the 
master equation Eq. (14) with the initial condition p n (0) = S nm . 

An explicit result for the noise in this case can be obtained by making 
further assumptions about the QDS and the coupling to the leads, see the 
following sections. For the general case, we only estimate AS. The current 
is of the order / ~ ew, with w some typical value of the cotunneling rate 
w nm , and thus 51 ~ ew. The time between switching from one dot-state to 
another due to cotunneling is approximately r ~ The correction AS" to 
the Poissonian noise can be estimated as AS" ~ 5I 2 t ~ e 2 w , which is of the 
same order as the Poissonian contribution el ~ e 2 w. Thus the correction 
to the Fano factor is of order unity. (Note however, that under certain 
conditions the Fano factor can diverge, see Sees. 5 and 6.) In contrast to this, 
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we find that for elastic cotunneling the off-diagonal rates vanish, w nm oc 
5 nm , and therefore 5p nn = and AS* = 0. Moreover, at zero temperature, 
either w nn (2, 1) or w nn (l, 2) must be zero (depending on the sign of the bias 
Ap). As a consequence, for elastic cotunneling we find Poissonian noise, 
F = S(0)/e\I\ = 1. 



5. Cotunneling through nearly degenerate states 

Suppose the QDS has nearly degenerate states with energies E n , and level 
spacing SE nm = E n — E m , which is much smaller than the average level 
spacing SE. In the regime, Ap,kBT,5E nm <C SE, the only allowed cotun- 
neling processes are the transitions between nearly degenerate states. The 
noise power is given by Eqs. (18) and (19), and below we calculate the 
correlation correction to the noise, AS. To proceed with our calculation we 
rewrite Eq. (14) for Sp(t) as a second-order differential equation in matrix 
form 

5p(t)=W 2 5p(t), Sp(0) = l-p, (20) 

where W is defined as W nm = w nm — 5 nm J2m' w m'n- We solve this equation 
by Fourier transformation, 

2W 

s ^> = "WOT (21) 

where we have used Wp = 0. We substitute Sp from this equation into 
Eq. (19) and write the result in a compact matrix form, 



A5H = -e 2 ^ 



(22) 



This equation gives the formal solution of the noise problem for nearly 
degenerate states. As an example we consider a two-level system. 

Using the detailed balance equation, W21P1 = we obtain for the 

stationary probabilities pi = ^12/(^12 + ^21), and p2 = w 21/(^12 + W2\)- 
From Eq. (17) we get 

j = c ^i2(w{ 1 +w I 21 ) + w 2 i(w I 22 + w[ 2 ) 

W 12 + U>21 

A straightforward calculation with the help of Eq. (21) gives for the cor- 
rection to the Poissonian noise 



A5(w) 



2e 2 {w I ll +w I 21 



U>22 ~ W 12) 



[w\2 + W2l)[u 2 + {w\2 + W2l) 2 \ 

x [w[ lWl 2W 2 i + w[ 2 w\ x -(!<-> 2)] . (24) 
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In particular, the zero frequency noise A5(0) diverges if the "off-diagonal" 
rates w nm vanish. This divergence has to be cut at to, or at the relaxation 
rate Wi n due to coupling to the bath (since w±2 in this case has to be 
replaced with w\2 + w- m ). The physical origin of the divergence is rather 
transparent: If the off-diagonal rates ^12,^21 are small, the QDS goes into 
an unstable state where it switches between states 1 and 2 with different 
currents in general [12]. The longer the QDS stays in the state 1 or 2 the 
larger the zero-frequency noise power is. However, if + w?2i = ^22 + w i2' 
then AS(lo) is suppressed to 0. For instance, for the QDS in the spin- 
degenerate state with an odd number of electrons AS(lj) = 0, since the two 
states I I) and | j) are physically equivalent. The other example of such a 
suppression of the correlation correction AS to noise is given by a multi- 
level QDS, 5E <C Ec, where the off-diagonal rates are small compared to 
the diagonal (elastic) rates [2]. Indeed, since the main contribution to the 
elastic rates comes from transitions through many virtual states, which do 
not participate in inelastic cotunneling, they do not depend on the initial 
conditions, = w^, and cancel in the numerator of Eq. (24), while they 
are still present in the current. Thus the correction AS /I vanishes in this 
case. Further below in this section we consider a few-level QDS, 5E ~ Eq, 
where AS / 0. 

To simplify further analysis we consider for a moment the case, where 
the singularity in the noise is most pronounced, namely, lo = and | <5E7i2 1 *C 
A/j,, k B T, so that w[ 2 = w^i, and W12 = W21. Then, from Eqs. (23) and (24) 
we obtain 

I = \(h+h), In = e (25) 

m=l,2 

A5(0) = (/l ~ 12)2 , (26) 

where I n is the current through the n-th level of the QDS. Thus in case 
\SE\2\ <C Afj,,k B T the following regimes have to be distinguished: (1) If 
k B T < Afi, then I n oc A/x, w\2 oc Afi, and thus both, the total current 
/ = e~ 1 GoA^i, and the total noise S = FGdA/j, are linear in the bias 
Afj, (here Gd is the conductance of the QDS). The total shot noise in this 
regime is super-Poissonian with the Fano factor F ~ I/iewu) S> 1. (2) In 
the regime A/x < ksT < F 1 / 2 A/j, the noise correction (26) arises because 
of the thermal switching the QDS between two states n = 1,2, where the 
currents are linear in the bias, I n ~ GoAfi/e. The rate of switching is 
W12 oc k B T, and thus AS ~ FG D Afi 2 /(k B T). Since k B T/Afi < F 1 / 2 , the 
noise correction AS is the dominant contribution to the noise, and thus the 
total noise S can be interpreted as being a thermal telegraph noise [13]. (3) 
Finally, in the regime i ?1 / 2 A / u < k B T the first term on the rhs of Eq. (18) 



cotunneling.tex; 1/02/2008; 23:08; p. 10 



SHOT NOISE OF COTUNNELING CURRENT 



11 




Figure 2. Left: Double-dot (DD) system containing two electrons and being weakly 
coupled to metallic leads 1, 2, each of which is at the chemical potential fi\, /j.2- The 
tunneling amplitudes between dots and leads are denoted by T. The tunneling (td) 
between the dots results in a singlet-triplet splitting J ~ t d /U with the singlet being 
a ground state. The tunneling path between dots and leads 1 and 2 forms a closed 
loop (shown by arrows) so that the Aharonov-Bohm phase 4> will be accumulated by an 
electron traversing the DD. Right: The Fano factor F = S(lo)/I, with the noise power 
S(u) given in Eqs. (18) and (37), is plotted as a function of the Aharonov-Bohm phase 
4> for the normalized bias v = Afj,/ J — 2 and for four different normalized frequencies 
Q = uj/[G(2Afi — J)] — 0.1, 0.25, 0.5, and 1. Inset: the same, but with fixed frequency 
Q — 0.1, where the bias v takes the values 1.5, 3, and oo. 

is the dominant contribution, and the total noise becomes an equilibrium 
Nyquist noise, S = 2Gb^bT ■ 

6. Noise of double-dot system: Two-particle Aharonov-Bohm effect 

We notice that for the noise power to be divergent the off-diagonal rates 
w\2 and u>2i have to vanish simultaneously. However, the matrix w nm is 
not symmetric since the off-diagonal rates depend on the bias in a different 
way. On the other hand, both rates contain the same matrix element of 
the cotunneling amplitude (Dj,Di>), see Eqs. (15) and (16). Although in 
general this matrix element is not small, it can vanish because of differ- 
ent symmetries of the two states. To illustrate this effect we consider the 
transport through a double-dot (DD) system (see Ref. [6] for details) as an 
example. Two leads are equally coupled to two dots in such a way that a 
closed loop is formed, and the dots are also connected, see Fig. 2. Thus, in 
a magnetic field the tunneling is described by the Hamiltonian Eq. (3) with 

D^YsTljcUis, 1,3 = 1,2, (27) 

T n = T 22 = T* 2 = T* x = e^/ 4 T , (28) 

where the last equation expresses the equal coupling of dots and leads and 
(j) is the Aharonov-Bohm phase. Each dot contains one electron, and weak 
tunneling td between the dots causes the exchange splitting [14] J ~ t d /U 
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(with U being the on-site repulsion) between one spin singlet and three 
triplets 

\ S ) = ^[4^1- 4/^0), 

\T ) = ^=[4 ] 4 l +4 l 4 ] ]\0), (29) 

|T + ) = d t 1T 4 T |0), |T_) = d t u 4 i |0). 

In the case of zero magnetic field, <p = 0, the tunneling Hamiltonian V 
is symmetric with respect to the exchange of electrons, 1 <-> 2. Thus the 
matrix element of the cotunneling transition between the singlet and three 
triplets (S\V(E — Ho)~ 1 V\Ti) , i = 0, ±, vanishes because these states have 
different orbital symmetries. A weak magnetic field breaks the symmetry, 
contributes to the off-diagonal rates, and thereby reduces noise. Next, we 
consider weak and strong cotunneling regimes. 



6.1. WEAK COTUNNELING 



In this regime, I/e <C wi n , according to the non-equilibrium FDT (see 
Sec. 3) the zero-frequency noise contains the same information as the av- 
erage current (the Fano factor F = 1). Therefore, we first concentrate on 
current. We focus on the regime, A/x S> J, where inelastic cotunneling [15] 
occurs with singlet and triplet contributions being different, and where 
we can neglect the dynamics generated by J compared to the one gener- 
ated by the bias ("slow spins"). Close to the sequential tunneling peak, 
A_ <C A + ~ U, we keep only the term D\D2 in the amplitude (7). After 
some calculations we obtain 

I = e^C^GA/i, (30) 

Is') + cos (p(4 s 'dis4sd2s>) , (31) 

s,s' 

where G = 7r(e^T 2 /A_) 2 is the conductance of a single dot in the cotun- 
neling regime [16], and we assumed Fermi liquid leads with the tunneling 
density of states v. Eq. (30) shows that the cotunneling current depends 
on the properties of the equilibrium state of the DD through the coherence 
factor C{ip) given in (31). The first term in C is the contribution from 
the topologically trivial tunneling path (phase-incoherent part) which runs 
from lead 1 through, say, dot 1 to lead 2 and back. The second term (phase- 
coherent part) in C results from an exchange process of electron 1 with 
electron 2 via the leads 1 and 2 such that a closed loop is formed enclosing 
an area A (see Fig. 2). Note that for singlet and triplets the initial and 
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final spin states are the same after such an exchange process. Thus, in the 
presence of a magnetic field B, an Aharonov-Bohm phase factor ip = ABe/h 
is acquired. 

Next, we evaluate C(<p) explicitly in the singlet-triplet basis (29) and 
discuss the applications to the physics of quantum entanglement (see the 
Ref. [6]). Note that only the singlet \S) and the triplet \Tq) are entangled 
EPR pairs while the remaining triplets are not (they factorize). Assuming 
that the DD is in one of these states we obtain the important result 

C(ip) = 2 T cos ip . (32) 

Thus, we see that the singlet (upper sign) and the triplets (lower sign) 
contribute with opposite sign to the phase-coherent part of the current. One 
has to distinguish, however, carefully the entangled from the non-entangled 
states. The phase-coherent part of the entangled states is a genuine two- 
particle effect, while the one of the product states cannot be distinguished 
from a phase-coherent single-particle effect. Indeed, this follows from the 
observation that the phase-coherent part in C factorizes for the product 
states T± while it does not so for S, Tq. Also, for states such as | fj,) the 
coherent part of C vanishes, showing that two different (and fixed) spin 
states cannot lead to a phase-coherent contribution since we know which 
electron goes which part of the loop. 

Finally, we present our results [6] for the high-frequency noise in the 
quantum range of frequancies, lo ~ A_ <C A+, and in the slow-spin regime 
Afj, 3> J. This range of frequancies is beyond the regime of the applicabil- 
ity of the non-equilibrium FDT, and therefore there is no simple relation 
between the average current and the noise (see the Sec. 3). After lengthy 
calculations using the perturbation expansion of (6) up to third order in V 
we obtain 



S(u) = (evrz/T 2 ) 2 [X u + X! 

ImXuj = ^ir _ w) _ e ^ 2 ~ ^ ' (33) 

HeX w = — sign(m — /X2 + a;) in | 



2-KUJ ^1^2 

1 m i^ 2 ~ u \ , at \ 1 1 Ml ^ 1 

v(u — /xij in I 1 + v(u — /j,2) in 



2-KL) 



M2 Ml 



(34) 



where Hi = A_(Z). Thus the real part of S(u>) is even in to, while the 
imaginary part is odd. A remarkable feature here is that the noise acquires 
an imaginary (i.e. odd-frequency) part for finite frequencies, in contrast 
to single-barrier junctions, where ImS(uj) always vanishes since we have 
5I\ = —512 for all times. In double-barrier junctions considered here we 
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find that at small enough bias Afi <C A_ = (A_(l) + A_(2)) /2, the odd 
part, lmS(uj), given in (33) exhibits two narrow peaks at w = ±/x, which 
in real time lead to slowly decaying oscillations, 

S odd (t) = Tr(evT 2 ) 2 C(v) Sm f^ t/2) sin(tA_). (35) 

These oscillations again depend on the phase-coherence factor C with the 
same properties as discussed before. These oscillations can be interpreted 
as a temporary build-up of a charge-imbalance on the DD during an un- 
certainty time ~ 1/A_, which results from cotunneling of electrons and an 
associated time delay between out- and ingoing currents. 



6.2. STRONG COTUNNELING 



The fact that in the perturbation V all spin indices are traced out helps 
us to map the four-level system to only two states \S) and \T) classified 
according to the orbital symmetry (since all triplets are antisymmetric in 
orbital space) . In Appendix B we derive the mapping to a two-level system 
and calculate the transition rates w nm (l',l) (n,m = 1 for a singlet and 
n, m = 2 for all triplets) using Eqs. (15) and (16) with the operators Di 
given by Eq. (27). Doing this we obtain the following result 

w nm (l, 1) = w nm (2, 2) = w nm (l, 2) = 0, 



( (l+cos0)A// (l-cos0)(A//+J) 
\ 3(l-cos0)(A/x-J) 3(l+cos0)A / u 



(36) 



which holds close to the sequential tunneling peak, A_ <C A + ~ U (but 
still A_ » J, A/i), and for A/x > J. We substitute this equation into 
the Eq. (24) and write the correction A5(w) to the Poissonian noise as 
a function of normalized bias v = A/i/ J and normalized frequency f2 = 
euj/[G(2An -J)} 

A5( W ) = ^ ^-f 1 ^- 1 '-^' 1 --^ . (37) 
(2f - 1) J [S2^ + (1 — cos^J^J 

From this equation it follows that the noise power has singularities as a 
function of uj for zero magnetic field, and it has singularities at (j) = 27rm 
(where m is integer) as a function of the magnetic field (see Fig. 2). We 
would like to emphasize that the noise is singular even if the exchange 
between the dots is weak, J <C A/Lt. In the case A/i < J the transition from 
the singlet to the triplet is forbidden by conservation of energy, W2i(2, 1) = 
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0, and we immediately obtain from Eq. (24) that AS(u>) = 0, i.e. the total 
noise is Poissonian (as it is always the case for elastic cotunneling) . In 
the case of large bias, A/x 3> J, two dots contribute independently to the 
current I = 2e~ 1 GA / u, and from Eq. (37) we obtain the Fano factor 

_ 3 cos 2 (?i>(l — cos d>) . . . 

8 iV + (1 — coscp)^ 

This Fano factor controls the transition to the telegraph noise and then to 
the equilibrium noise at high temperature, as described above. We notice 
that if the coupling of the dots to the leads is not equal, then w nm (l, I) / 
serves as a cut-off of the singularity in AS (to). 

Finally, we remark that the Fano factor is a periodic function of the 
phase 4> (see Fig. 2); this is nothing but an Aharonov-Bohm effect in the 
noise of the cotunneling transport through the DD. However, in contrast 
to the Aharonov-Bohm effect in the cotunneling current through the DD 
which has been discussed earlier in the Sec. 6.1, the noise effect does not 
allow us to probe the ground state of the DD, since the DD is already in a 
mixture of the singlet and three triplet states. 

7. Cotunneling through continuum of single-electron states 

We consider now the transport through a multi-level QDS with 5E <C Eq- 
In the low bias regime, Afi <C (5E Eq) 1 ^ 2 , the elastic cotunneling dom- 
inates transport [2], and according to the results of Sec. 4 the noise is 
Poissonian. Here we consider the opposite regime of inelastic cotunneling, 
Afi S> (5E Ec) 1 / 2 ■ Since a large number M of levels participate in trans- 
port, we can neglect the correlations which we have studied in Sees. 5 and 
6, since they become a 1/M-effect. Instead, we concentrate on the heating 
effect, which is not relevant for the 2- level system considered before. The 
condition for strong cotunneling has to be rewritten in a single-particle 
form, Ti n 3> t c , where T m is the single-particle energy relaxation time on 
the QDS due to the coupling to the environment, and r c is the time of the 
cotunneling transition, which can be estimated as r c ~ evoAfi/I (where 
vn is the density of QDS states). Since the energy relaxation rate on the 
QDS is small, the multiple cotunneling transitions can cause high energy 
excitations on the dot, and this leads to a nonvanishing backward tunneling, 
w nm (l,2) / 0. In the absence of correlations between cotunneling events, 
Eqs. (17) and (18) can be rewritten in terms of forward and backward 
tunneling currents I + and i_ , 

/ = /+-/_, S = e(I+ + I-), (39) 
1+ = e^u/ nm (2, l)p m , 7_ = w nm (l,2)p m , (40) 
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where the transition rates are given by (15). 

It is convenient to rewrite the currents I± in a single-particle basis. 
To do so we substitute the rates Eq. (15) into Eq. (40) and neglect the 
dependence of the tunneling amplitudes Eq. (3) on the quantum numbers 
k and p, Tik p = Tj, which is a reasonable assumption for QDS with a large 
number of electrons. Then we define the distribution function on the QDS 
as 

f{e) = v D 1 Y,Ke-e P )TTpdld p (41) 

v 

and replace the summation over p with an integration over e. Doing this 
we obtain the following expressions for T = 



GiG 2 / 1 1 



2 



c± = aT? // <fe<fe ' ( £ - £>± A ")/( S )I 1 - ( 43 > 

where Gi 2 = 7re 2 i/i/£)|Ti 2I 2 are the tunneling conductances of the barriers 
1 and 2, and where we have introduced the function 0(e) = e6(e) with 
6(e) being the step-function. In particular, using the property 0(e + A/i) — 
0(e — A/x) = e + A[i and fixing 

J de[f{e) - 9(-e)} = 0, (44) 

(since I± given by Eq. (42) and Eq. (43) do not depend on the shift e — > 
e + const) we arrive at the following general expression for the cotunneling 
current 



2 



A = 1 + 12T/(A^) 2 , (46) 
dee[f(e) - 9{-e)} > 0, (47) 



where the value has the physical meaning of the energy acquired by 
the QDS due to the cotunneling current through it. 

We have deliberately introduced the functions C± in the Eq. (42) to 
emphasize the fact that if the distribution /(e) scales with the bias An (i.e. 
/ is a function of e/A/x), then C± become dimensionless universal numbers. 
Thus both, the prefactor A (given by Eq. (46)) in the cotunneling current, 
and the Fano factor, 
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take their universal values, which do not depend on the bias Afi. We con- 
sider now such universal regimes. The first example is the case of weak 
cotunneling, T Ui <C r c , when the QDS is in its ground state, /(e) = 9(—s), 
and the thermal energy of the QDS vanishes, T = 0. Then A = 1, and 
Eq. (45) reproduces the results of Ref. [2]. As we have already mentioned, 
the backward current vanishes, /_ = 0, and the Fano factor acquires its 
full Poissonian value F = 1, in agreement with our nonequilibrium FDT 
proven in Sec. 3. In the limit of strong cotunneling, T m S> r c , the energy 
relaxation on the QDS can be neglected. Depending on the electron-electron 
scattering time r ee two cases have to be distinguished: The regime of cold 
electrons r ee 3> r c and regime of hot electrons r ee <C r c on the QDS. Below 
we discuss both regimes in detail and demonstrate their universality. 

7.1. COLD ELECTRONS 

In this regime the electron-electron scattering on the QDS can be neglected 
and the distribution /(e) has to be found from the master equation Eq. (14). 
We multiply this equation by Yl p $( £ — e v) (n\d P d p \n) , sum over n and 
use the tunneling rates from Eq. (15). Doing this we obtain the standard 
stationary kinetic equation which can be written in the following form 



where A = (G\ + G\) / {2G1G2) > 1 arises from the equilibration rates 
w mn (l,l). (We assume that if the limits of the integration over energy e 
are not specified, then the integral goes from —00 to +00.) From the form 
of this equation we immediately conclude that its solution is a function 
of e/A/x, and thus the cold electron regime is universal as defined in the 
previous section. It is easy to check that the detailed balance does not hold, 
and in addition a{e) 7^ a(—e). Thus we face a difficult problem of solving 
Eq. (49) in its full nonlinear form. Fortunately, there is a way to avoid this 
problem and to reduce the equation to a linear form which we show next. 
We group all nonlinear terms on the rhs of Eq. (49): J de'a (e'— e)/(e') = 



h(e)f(e), where h(e) = J de' {a(e' - e)f(e') + a(e - e')[l - f(e')}} . The 



trick is to rewrite the function h{e) in terms of known functions. For doing 
this we split the integral in h{e) into two integrals over e' > and e' < 0, 
and then use Eq. (44) and the property of the kernel a(e)—a(—e) = 2(1+A)e 
to regroup terms in such a way that h(e) does not contain /(e) explicitly. 




(49) 



(50) 



± 
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Taking into account Eq. (47) we arrive at the following linear integral 
equation 

/ de'a{e' - e)f(e') = [(1 + X)(e 2 + 2T) + (A^) 2 ]/(e), (51) 

where the parameter T is the only signature of the nonlinearity of Eq. (49) . 

Since Eq. (51) represents an eigenvalue problem for a linear operator, 
it can in general have more than one solution. However, there is only one 
physical solution, which satisfies the conditions 

0</(e)<l, /(-oo) = l, /(+oo) = 0. (52) 

Indeed, using a standard procedure one can show that two solutions of the 
integral equation (51), f± and /2, corresponding to different parameters 
Yi / ^2 should be orthogonal, J cfe/i(e)/2(— e) = 0. This contradicts the 
conditions Eq. (52). The solution is also unique for the same T, i.e. it is 
not degenerate (for a proof, see the Ref. [5]). From Eq. (49) and conditions 
Eq. (52) it follows that if /(e) is a solution then 1 — /(— e) also satisfies 
Eqs. (49) and (52). Since the solution is unique, it has to have the symmetry 
f(e) = l-f(-e). 

2.2 
2.1 
2.0 
< 1.9 
1.8 
1.7 
1.6 

123456789 10 

G 1 /G 2 

Figure 3. The prefactor A in the expression (45) for the cotunneling current characterizes 
a universal cotunneling transport in the regime of weak cotunneling, r in <C t c , (A = 1, 
see Ref. [2]), and in the regime of strong cotunneling, Tin 3> t c (A > 1). Here A is plotted 
as a function of G1/G2 (same as a function of G2/G1) for the strong cotunneling, for the 
cold-electron case, r ee ^> r c (solid line) and for the hot-electron case, r ee <JC r c (dotted 
line). Gi,2 are the tunneling conductances of a junctions connecting leads 1 and 2 with 
the QDS. 

We solve Eqs. (51) and (52) numerically and use Eqs. (43) and (48) to 
find that the Fano factor is very close to 1 (it does not exceed the value 



cold electrons 
hot electrons 
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F 1.006). Next we use Eqs. (46) and (47) to calculate the prefactor A and 
plot the result as a function of the ratio of tunneling conductances, G1/G2, 
(Fig. 3, solid line). For equal coupling to the leads, G\ = G2, the prefactor 
A takes its maximum value 2.173, and thus the cotunneling current is 
approximately twice as large compared to its value for the case of weak 
cotunneling, T m <C r c . A slowly decreases with increasing asymmetry of 
coupling and tends to its minimum value A = 1 for the strongly asymmetric 
coupling case G1/G2 or G2/G1 3> 1. 

7.2. HOT ELECTRONS 

In the regime of hot electrons, r ee <C t c , the distribution is given by the 
equilibrium Fermi function /f(^) = [1 + exp(e/fesT e )] _1 , while the electron 
temperature T e has to be found self-consistently from the kinetic equation. 
Eq. (49) has to be modified to take into account electron-electron interac- 
tions. This can be done by adding the electron collision integral I ee {£) to 
the rhs of (49) . Since the form of the distribution is known we need only the 
energy balance equation, which can be derived by multiplying the modified 
equation (49) by e and integrating it over e. The contribution from the 
collision integral I ee (e) vanishes, because the electron-electron scattering 
conserves the energy of the system. Using the symmetry = 1 — /f( — e) 

we arrive at the following equation 

J J dede'f F (e')[l - f F (e)]a(e' - e)e = 0. (53) 

Next we regroup the terms in this equation such that it contains only 
integrals of the form f °° defF{e){- • •)• This allows us to get rid of nonlinear 
terms, and we arrive at the following equation, 

J dee 3 [f F (e) - 9(-e)} + 3T 2 = , (54) 

which holds also for the regime of cold electrons. Finally, we calculate the 
integral in Eq. (54) and express the result in terms of the dimensionless 
parameter a = A^//csT e , 

a = vr[8(l + A)/5] 1/4 . (55) 

Thus, since the distribution again depends on the ratio e/Afj,, the hot 
electron regime is also universal. 

The next step is to substitute the Fermi distribution function with the 
temperature given by Eq. (55) into Eq. (43). We calculate the integrals and 
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arrive at the closed analytical expressions for the values of interest, 



where again A = {G\ + G|) /2G1G2 > 1. It turns out that similar to the 
case of cold electrons, Sec. 7.1, the Fano factor for hot electrons is very 
close to 1 (namely, it does not exceed the value F f» 1.007). Therefore, we 
do not expect that the super-Poissonian noise considered in this section 
(i.e. the one which is due to heating of a large QDS caused by inelastic 
cotunneling through it) will be easy to observe in experiments. On the 
other hand, the transport-induced heating of a large QDS can be observed 
in the cotunneling current through the prefactor A, which according to 
Eq. (56) takes its maximum value A = 1 + \/5/4 ~ 2.118 for G± = G2 and 
slowly reaches its minimum value 1 with increasing (or decreasing) the ratio 
G1/G2 (see Fig. 3, dotted line). Surprisingly, the two curves of A vs G1/G2 
for the cold- and hot-electron regimes lie very close, which means that the 
effect of the electron-electron scattering on the cotunneling transport is 
rather weak. 

8. Conclusions 

Here we give a short summary of our results. In Sec. 3, we have derived the 
non-equilibrium FDT, i.e. the universal relation (12) between the current 
and the noise, for QDS in the weak cotunneling regime. Taking the limit 
T, lo — ► 0, we show that the noise is Poissonian, i.e. F = 1. 

In Sec. 4, we present the results of the microscopic theory of strong 
cotunneling, Ref. [5]: The master equation, Eq. (14), the average current, 
Eq. (17), and the current correlators, Eqs. (18) and (19), for a QDS system 
coupled to leads in the strong cotunneling regime w- m <C I/e at small 
frequencies, uj <C A±. In contrast to sequential tunneling, where shot noise 
is either Poissonian (F = 1) or suppressed due to charge conservation 
(F < 1), we find that the noise in the inelastic cotunneling regime can be 
super-Poissonian (F > 1), with a correction being as large as the Poissonian 
noise itself. In the regime of elastic cotunneling F = 1. 

While the amount of super-Poissonian noise is merely estimated at the 
end of Sec. 4, the noise of the cotunneling current is calculated for the 
special case of a QDS with nearly degenerate states, i.e. SE nm <C 5E, in 
Sec. 5, where we apply our results from Sec. 4. The general solution Eq. (22) 
is further analyzed for two nearly degenerate levels, with the result Eq. (24). 




(56) 



(57) 



cotunneling.tex; 1/02/2008; 23:08; p. 20 



SHOT NOISE OF COTUNNELING CURRENT 



21 



More information is gained in the specific case of a DD coupled to leads 
considered in Sec. 6, where we determine the average current Eqs. (30-31) 
and noise Eqs. (33-34) in the weak cotunneling regime and the correlation 
correction to noise Eq. (37) in the strong cotunneling regime as a function 
of frequency, bias, and the Aharonov-Bohm phase threading the tunneling 
loop, finding signatures of the Aharonov-Bohm effect and of the quantum 
entanglement. 

Finally, in Sec. 7, another important situation is studied in detail, the 
cotunneling through a QDS with a continuous energy spectrum, 5E <C 
A// <C Eq- Here, the correlation between tunneling events plays a minor 
role as a source of super-Poissonian noise, which is now caused by heating 
effects opening the possibility for tunneling events in the reverse direction 
and thus to an enhanced noise power. In Eq. (48), we express the Fano 
factor F in the continuum case in terms of the dimensionless numbers C±, 
defined in Eq. (43), which depend on the electronic distribution function 
/(e) in the QDS (in this regime, a description on the single-electron level is 
appropriate). The current Eq. (45) is expressed in terms of the prefactor A, 
Eq. (46). Both F and A are then calculated for different regimes. For weak 
cotunneling, we immediately find F = 1, as anticipated earlier, while for 
strong cotunneling we distinguish the two regimes of cold (r ee 3> r c ) and 
hot (r ee <C t c ) electrons. For both regimes we find that the Fano factor is 
very close to one, while A is given in Fig. 3. 
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Appendix A 

In this Appendix we present the derivation of Eqs. 7 and 8. In order to 
simplify the intermediate steps, we use the notation 0(t) = f^dt'Otf) 
for any operator O, and O(0) = O. We notice that, if an operator O is a 
linear function of operators Di and D, , then O(oo) = (see the discussion 
in Sec. 3). Next, the currents can be represented as the difference and the 
sum of I\ and 1%, 



where X = D2 + D\, and Y = D\ + D2. While for the perturbation we have 



Id = (/ 2 - A)/2 = iepft - X)/2, 
Is = (ii + / 2 )/2 = ie(Y* - Y)/2 , 



(58) 
(59) 



v = x + x ] = y + y f . 



(60) 
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First we concentrate on the derivation of Eq. (7) and redefine the average 
current Eq. (5) as I = Id (which gives the same result anyway, because the 
average number of electrons on the QDS does not change I s = 0). 

To proceed with our derivation, we make use of Eq. (6) and expand the 
current up to fourth order in T^ p : 

o t o 

I = i J dt j dt'{i d V(t)V(t')V{t')) -i J dt{VI d V{t)V{t)) + c.c. (61) 

— oo — oo — oo 

Next, we use the cyclic property of trace to shift the time dependence to 
Id- Then we complete the integral over time t and use I d (oo) = 0. This 
procedure allows us to combine first and second term in Eq. (61), 

o 

I = -i J dt{[I d V + vi d ]V(t)V(t))+ c.c. (62) 

— oo 

Now, using Eqs. (58) and (60) we replace operators in Eq. (62) with X 
and Xt in two steps: I = e f ^ dt{[X*X* - XX]V{t)V{t)) + c.c, where 
some terms cancel exactly. Then we work with V(t)V(t) and notice that 
some terms cancel, because they are linear in cik and c\ k . Thus we obtain 
I = e f^dtdX^X^ - XX][X^(t)X^(t) + X{t)X(t)]) + c.c. Two terms 
XXXX and X^X^X^X^ describe tunneling of two electrons from the same 
lead, and therefore they do not contribute to the normal current. We then 
combine all other terms to extend the integral to +oo, 

oo 

I = e J dt(X\t)X\t)XX - XXX\t)X\t)) (63) 

— oo 

Finally, we use dtX{t)X{t) = - dtX(t)X(t) (since X{po) = 0) to 

get Eq. (7) with A = XX. Here, again, we drop terms -Di-D} and D2D2 
responsible for tunneling of two electrons from the same lead, and obtain 
A as in Eq. (7). 

Next, we derive Eq. (8) for the noise power. At small frequencies u <C 
A± fluctuations of I s are suppressed because of charge conservation (see 
below), and we can replace I2 in the correlator Eq. (5) with I d . We expand 
S(lo) up to fourth order in T^ p , use f^™ dt Id{t)e ±lu)t = 0, and repeat the 
steps leading to Eq. (62). Doing this we obtain, 

00 

S(u) = - J dtcos(Ljt)([V(t),i d (t)][V,i d }). (64) 
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Then, we replace V and Id with X and X< . We again keep only terms 
relevant for cotunneling, and in addition we neglect terms of order u/ A± 
(applying same arguments as before, see Eq. (65)). We then arrive at Eq. (8) 
with the operator A given by Eq. (7). 

Finally, in order to show that fluctuations of I s are suppressed, we re- 
place Id in Eq. (64) with I s , and then use the operators Y and instead of 
X and X^ . In contrast to Eq. (63) terms such as Y^Y^YY do not contribute, 
because they contain integrals of the form J^^dtcos(wt)Di(t)Di'(t) = 0. 
The only nonzero contribution can be written as 

oo 

S ss (u) = ^f I dtco S (ut)([Y\t),Y(t)][Y\Y]), (65) 

— oo 

where we have used integration by parts and the property Y(oo) = 0. 
Compared to Eq. (8) this expression contains an additional integration 
over t, and thereby it is of order (u>/A±) 2 . 



Appendix B 



In this Appendix we calculate the transition rates Eq. (15) for a DD coupled 
to leads with the coupling described by Eqs. (27) and (28) and show that 
the four-level system in the singlet-triplet basis Eq. (29) can be mapped 
to a two-level system. For the moment we assume that the indices n and 
m enumerate the singlet-triplet basis, n,m = S,Tq,T + ,T-. Close to the 
sequential tunneling peak, A_ <C A+, we keep only terms of the form 
DjDi>. Calculating the trace over the leads explicitly, we obtain at T = 0, 



2A 1 



^2 T lj T lj' T l*j' T l'j M nm(j,f) , 



3,3 



M nm (j,j') = ^(n|d^d s 'j|m)(m|4, j .,d 

s,s' 



r n) 



(66) 
(67) 



with Q(e) = eO(e), and 5E nm = 0, ±J, and we have assumed tj, <C A_. 

Since the quantum dots are the same we get M nm (l, 1) = M nm (2, 2), 
and M nm (l,2) = M nm (2, 1). We calculate these matrix elements in the 
singlet-triplet basis explicitly, 

/ 1 1 
1 1 



M(l,l) 



1 1 



1\ 
1 





(68) 



\1 10 2/ 
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/ 


1 


-1 


-1 


-1\ 


1 




-1 


1 


1 


1 


2 




-1 


1 


2 







V 


-1 


1 





2 / 



(69) 



Assuming now equal coupling of the form Eq. (28) we find that for I = I' 
the matrix elements of the singlet-triplet transition vanish (as we have 
expected, see Sec. 5). On the other hand the triplets are degenerate, i.e. 
5E nm = in the triplet sector. Then from Eq. (66) it follows that w nm (l, I) = 
0. Next, we have 0(^2 — Mi — 5E nrn ) = 0, since for nearly degenerate states 
we assume A/i > |<5^nm|; and thus w nm (1,2) =0. Finally, for w nm = w nm = 
(2, 1) we obtain, 

2 

H'.s-.s - - ! '— - ) A//U • c(.su). (70) 
w st = 77 ( ) (A// + J)(l — cos <fr), (71) 
(r;/ ,, _ _ ( __ ) (A/x- J)(l-cos0), (72) 



WTT 




(73) 



Next we prove the mapping to a two-level system. First we notice that 
because the matrix wtt is symmetric, the detailed balance equation for the 
stationary state gives p n /Pm = w rnn /w nrn = 1, n, m G T. Thus we can set 
Pn — ► P2/3, for n e T. The specific form of the transition matrix Eqs. (70- 
73) helps us to complete the mapping by setting (1/3) ^^=2^1"* ~^ w i2, 
Y.n=2 w m -»■ w 2 \, and (1/3) ^ nm=2 w nm -»■ ^22, so that we get the 
new transition matrix Eq. (36), while the stationary master equation for 
the new two-level density matrix does not change its form. If in addi- 
tion we set (l/3)X)m=2<fy>i™(*) -> 6pi 2 (t), En=2<W(*) $P2i(t), and 
(1/3) m=2 $Pnm(t) — ► 8p22(t), then the master equation Eq. (14) for 
fipnm(t) and the initial condition (5p nm (0) = (5 nm — p n do not change either. 
Finally, one can see that under this mapping Eq. (19) for the correction 
to the noise power AS(lv) remains unchanged. Thus we have accomplished 
the mapping of our singlet-triplet system to the two-level system with the 
new transition matrix given by Eq. (36). 
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